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ABSTRACT 

We have cross-correlated the WMAP data with several surveys of extragalactic sources 
and find evidence for temperature decrements associated with galaxy clusters and 
groups detected in the APM Galaxy Survey survey and the ACO catalogue. We 
interpret this as evidence for the thermal Sunyaev-Zeldovich (SZ) effect from the 
clusters. Most interestingly, the signal may extend to sal deg (« 5/i _1 Mpc) around 
both groups and clusters and we suggest that this may be due to hot 'superclus- 
ter' gas. We have further cross-correlated the WMAP data with clusters identified 
f*^ ' in the 2MASS galaxy catalogue and also find evidence for temperature decrements 

OO ! there. From the APM group data we estimate the mean Compton parameter as 

y(z < 0.2) = 7 ± 3.8 x 10 -7 . We have further estimated the gas mass associated 
^sO ' with the galaxy group and cluster haloes. Assuming temperatures of 5 keV for ACO 

, clusters and 1 keV for APM groups and clusters, we derive average gas masses of 

M{r < 1.75/i _1 Mpc) r* 3 x 10 13 h~ 2 M Q for both, the assumed gas temperature and 
SZ central decrement differences approximately cancelling. Using the space density of 
APM groups we then estimate n gas « 0.03k- 1 (1 keV/kT)(9 max /20') - 75 . For an SZ 
P^i 1 extent of 6 max = 20', kT = 1 keV and h = 0.7, this value of Q, gas w 0.04 is consistent 

with the standard value of Q barvon = 0.044 but if the indications we have found for a 
^ ■ more extended SZ effect out to 9 max w 60' are confirmed, then higher values of Q,Q 9as 

^ | will be implied. Finally, the contribution to the WMAP temperature power spectrum 

C$ from the extended SZ effect around the z < 0.2 APM+ACO groups and clusters is 

1-2 orders of magnitude lower than the I = 220 first acoustic peak. But if a similar 
SZ effect arises from more distant clusters then this contribution could increase by a 
. factor > 10 and then could seriously affect the WMAP cosmological fits. 

S-H 

' Key words: cosmic microwave background - galaxies: clusters 



1 INTRODUCTION 

The Wilkinson Microwave Anisotropy Probe (WMAP) cos- 
mic microwave background (CMB) anisotropy experiment 
has published its first year data and has already provided 
remarkable results. It has confirmed that the first peak in 
the power-spectrum occurs at / = 220 ± 10 and has pro- 
vided an excellent f urther detection of the second peak 
jHinshaw et alJ [2003) . In the first instance, these new re- 
sults seem to pro vide further suppo rt for the standard 
ACDM cosmology dSpergel et al.ll2003h . But the main new 
result from WMAP is the detection of polarisation at large 
scales which can only ha ve come from an e poch of reion- 
isation at 10 < z < 20 iKoeut et a!ll2003fi . This results 
in a significant, («30 per cent), reduction of the acoustic 
peak heights due to Thomson scattering and shows that 
the CMB anisotropies are seriously affected by low red- 
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shift galaxy formation physics. Furthermore, WMAP also 
finds a low quadrupo le in the temperature power-spectrum 
Jffinshaw et alj l2003) and this is also unexpected if the In- 
tegrated Sachs- Wolfe (ISW) effect caused by relatively re- 
cent domination of the cosmological constant is present. 
Therefore, although WMAP data supports the standard 
model, there are significant problems and already strong ev- 
idence that the temperature spectrum is affected by cosmic 
foregrounds such as the re-ionised intergalactic medium at 
10 < z < 20. 

Here we investigate the possibility that other low red- 
shift processes have filtered the WMAP temperature spec- 
trum. In particular, we shall search for Sunyaev-Zeldovich 
(SZ), inverse Compton-scattering of microwave background 
photons by hot gas in clusters. Various authors have made 
model dependent predictions for contamination of the CMB 
data from the SZ effect and usu ally concluded that t he 
conta minatin g effects were s mall jRefregier et alJ l2000allbl 
iKomatsu fc Kitavamalll999l) . The new WMAP data gives 
a first real opportunity to make empirical checks of the level 
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of contamination by direct cross-correlation of the high reso- 
lution 94GHz W band with galaxy cluster data. The WMAP 
team themselves have presented evidence for both SZ and ra- 
dio so urce contamination in the WMAP data. [Bennett et alJ 
( 2003) list 208 point sources detected at more than 5a in 
the WMAP data and identified them as radio galaxies and 
quasars. They also find significant WMAP W band SZ detec- 
tions of the brightest X-ray clusters such as Coma and also 
found a 2.5a detection of the SZ effect in the XBACS sam - 
ple of 242 X-ray bright Abell clusters (lEbeling et alJll99rj) . 
However, they only looked for SZ decrements on the scale 
of the WMAP beam and did not explore any larger scales. 

While this paper was in preparation, several other re - 
late d papers have app eared. iBouehn fc Crittenden! l|2003l) 
and iNolta et alJ i2003T) have claimed evidence for the ISW 
effect in the WM AP data from cross- c orrela tion with the 
NVSS catalogue. iDieeo. Silk fc Sliwal (|2003h have cross- 
correlated X-ray dat a and the WMAP data. Similarly 
iFosalba fc Gaztanaeal {2(303) have cross-correlated APM 
galaxies and WMAP data; they find a marginal detection of 
the ISW effect at 5-10 deg scales and suggest that the lack of 
a detection at sm aller scales may be due to cancellation with 
an S Z effect. iHernandez-Monteaeudo fc Rubino-Martml 
( 2003) have obtained upper limits on diffuse SZ emis- 
sion from superclusters from failing to find any cross- 
correlation with the Abell-Corwin-Olowin (ACO) and other 
cluster catalogues but the y did find significant corre l ations 
from individual clusters. iGiommi fc Colafrancescol (120031) 
have claimed to detect significa nt blazar contamination 
in WM AP and Boomerang data. Afs hordi. Loh fc Strauss! 
( 2003) have claimed the detection of SZ, ISW and point 
sources in a power spectrum analysis of WMAP data and 
the 2MASS galaxy catalogue. 



2 DATA 
2.1 WMAP 

The W MAP data has been published by iBennett et alJ 
(2003) in HEALPix format. Here we shall principally use 
the WMAP W band data at 94 GHz because of its rela- 
tively high resolution compared to the other bands. The 
FWHM of the 94 GHz W beam is 12.' 6 compared to 19/8 at 
V (61 GHz), 29/4 at Q (41 GHz), 37/2 at Ka (33 GHz) and 
49/2 at K (23 GHz). Althoug h, none of the beams is exactly 
Gaussian (see Fig. 2 of lPage et alJ (|2003r0 we have found 
that Gaussians of the above FWHM are good approximate 
fits to cross-correlation results between faint radio point 
sources and WMAP data. We shall also use the Internal 
Linear Combination (ILC) map which has 1 deg resolution 
(Benn ett et alj|2003f) . Also the W band has most sensitivity 
to the SZ effect via its high resolution, although perhaps has 
higher noise than the V band. The HEALPix WMAP data 
we shall use has equal area pixels of 49 arcmi n 2 . Where nec- 
essary , we shall use the KpO WMAP mask of lBennett et alJ 
( 2003) which mainly masks Galactic contamination but usu- 
ally its effect is small because we shall generally be working 
at Galactic latitudes, \b\ >40 deg. The maps all use thermo- 
dynamic temperature and the cosmological dipole has al- 
ready been subtracted from the data by the WMAP team. 



2.2 Galaxy Cluster Catalogues 

We shall be using thr ee galaxy cluster catalogues. The first is 
the ACO catalogue of lAbell. Corwin fc OlowirJ Il989ft which 
lists clusters with 30 or more members in a 1.5k- 1 Mpc ra- 
dius within 2 mag of the 3rd brightest cluster member. It 
assigns richness class (0 ^ R ^ 5) with Coma classed as 
R = 3. The Northern catalogue with b > 40 deg lists 2489 
clusters with R ^ the Southern catalogue with b < 40 
deg lists 1346 such clusters. The sky density in the North is 
therefore 0.52 deg" 2 and in the South it is 0.28 deg" 2 . The 
sky density of R ^ 2 clusters is 0.063 deg -2 with an average 
redshift of z = 0.15. 

We shall also use galaxy group and clu ster catal ogues 
derive d from the APM Galaxy Survey of iMaddox et alJ 
(1990) which covers the whole area with 8 < —2.5 deg and 
b < —40 deg. These wer e identified u s ing th e same 'friends- 
of-friends' algorithm as iMvers et al.l {2,003) and references 
therein. Circles around each APM galaxy with B < 20.5 are 
'grown' until the overdensity, (3, falls to j3 = 8 and those 
galaxies whose circles overlap are called groups. Minimum 
memberships, m, of m 7 and m ^ 15 were used which de- 
fine minimum group effective 'radii' of 1/2 and 1/7, since the 
APM galaxy surface density is N » 750 deg" 2 at B < 20.5. 
We assume an average redshift of z = 0.1 for both APM 
samples. The sky density of groups and clusters is 3.5 deg -2 
at m 7 and 0.35 deg -2 at m ^ 15. Even at m ^ 15 there 
are differences between the ACO and APM catalogues; an 
R — Abell cluster at z — 0.15 may contain only a galaxy 
sky density of 260 deg" 2 within its 11/5 Abell radius com- 
pared to a minimum sky density of 5250 deg -2 for galaxies 
within the APM groups. Therefore the R < 1 ACO 'clusters' 
may accommodate much lower density galaxy associations 
than even the m ^ 7 APM groups which are guaranteed to 
sample higher density regions, albeit over smaller areas. 

The third cluster catalogue is derived from the final data 
release of the 2MA SS Extended Source Catalogue (XSC) 
dJarrett et alJ l200d) to a limit of K s < 13.7. ^-selected 
galaxy samples are dominated by early-type galaxies be- 
cause of their red colours and early-type galaxies are the 
most common galaxy-type found in rich galaxy clusters. 
Therefore the 2MASS survey provides an excellent tracer 
of the high density parts of the Universe out to z < 0.15 
and so provides a further test for the existence of the SZ 
effect. Using the above 2-D friends-of-friends algorithm, we 
have detected 500 groups and clusters with m 35 mem- 
bers at the density contrast j3 — 8 in the |6| ^ 10 deg area. 
The 2MASS groups have average redshift, z w 0.06. 



3 CROSS-CORRELATION TECHNIQUE 

We search for correlated signals between the 5 (+ ILC) 
WMAP bands and the above datasets using a simple cross- 
correlation technique. The WMAP data are presented in the 
form of temperature differences with respect to the global 
mean and we form the correlation function by calculating 
the average AT in annuli around each cluster as a function 
of angular distance, 9, between the cluster and pixel centres. 

At large 8, where the average background may be ex- 
pected to be sampled, the average in particular areas of the 
sky may not return to zero, probably because of Galactic 
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contamination. We therefore always calculate the simple av- 
erage of all temperature differences in the sub-area of sky 
used and plot this on the cross-correlation result as the ef- 
fective zero level (solid lines). 

The errors are calculated by repeated Monte Carlo ran- 
dom realisations of the cluster centres across the WMAP 
area used. The realisations always have the same number of 
points as the parent sample and 100 realisations were done 
in each case. The effect of clustering of clusters is included 
in the realisations; the clustering amplitude is matched to 
that of the cluster sample under consideration. In the cross- 
correlation functions, there is significant correlation between 
the results at different angles, 9, so we have also calculated 
Monte-Carlo errors in the integrated bins used to quote over- 
all significances. The Monte Carlo realisations also quite ac- 
curately return results that are close to the overall average, 
indicating that our cross-correlation technique is free of any 
systematics. We note that our error analysis will miss vari- 
ance from CMB modes of order the cluster catalogue size. 
Since this size is generally > 50 deg, it is not expected that 
error estimates on angular scales 6 < 1 deg will be seriously 
affected, although on larger scales our quoted errors may 
increasingly be underestimates. 



4 CROSS-CORRELATION RESULTS 

In Fig. la,b we first present the results from cross-correlating 
our APM groups and cluster centres with the WMAP W 
band data. We consider the two APM cluster samples with 
m ^ 15 and m ^ 7 members. Both samples show indications 
of anti-correlation with respect to the 94 GHz data. Relative 
to the overall mean, the integrated significances for the m ^ 
15 sample range from 1.7a at 9 < 10' to 1.0a at 9 < 30'. 
The significances for the m ^ 7 sample range from 2.0a 
at 9 < 10' to 1.8cr at 9 < 30'. In the case of the m ^ 15 
clusters, the size of the anti-correlation is AT ~ — O.OlmK 
on WMAP pixel scales. It is also interesting that the size 
of the signal only goes down slightly to AT ~ —0.008 mK, 
at small scales for the m ^ 7 sample and here the anti- 
correlation appears to extend beyond the beam-size, only 
reducing to 1.5a at 9 ~ 60'. The anti-correlation signal seen 
at small 9 10') is consistent with what is expected at this 
freque ncy for the temperatu re decrement caused by the SZ 
effect llRefreeier et a l. 2000a). The surprise is that it can be 
seen out to ml deg scales and that it persists in groups with 
a sky density of 3.5 deg -2 , a factor «30x higher than that 
of the Abell clusters. 

We now turn to the ACO catalogue to see if there is 
any confirmation of this tentative detection of extended anti- 
correlation around the APM groups and clusters. When the 
full ACO catalogue at |6| > 40 deg including clusters of all 
richnesses, was used only an insignificant anti-correlation 
was found. The catalogue was then cut back to R ^ 2 
and the result changed dramatically and consistently in 
both Hemispheres. Based on the remaining 229 clusters with 
b > 40 deg and the remaining 377 clusters with b < —40 deg, 
significant anti-correlation is again seen out to scales of 
9 « ldeg or mLSlr 1 Mpc (see Fig. lc). The Monte-Carlo 
errors indicate that the effect is significant at the ~ 2. la- 
level for 9 < 10' and at the ~ 2.2a level for 9 < 1 deg in the 
N+S sample. Our interpretation is that the R = 0, 1 clus- 
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Figure 1. (a) Cross-correlation of WMAP 94 GHz W band data 
with m ^ 15 APM groups and clusters, (b) Cross-correlation 
of WMAP 94 GHz W band data with m ^ 7 APM groups and 
clusters. The dashed line is an isothermal SZ model with AT(0) = 
0.015 mK, 8 C = 2/3, (3 = 0.75 convolved with a Gaussian beam 
profile with a = 5. '25. (c) Cross-correlation of WMAP 94 GHz 
W band data with ACO R > 2 clusters for combined ACO |6| > 
40 deg N+S samples. The dashed, dotted and dot-dash lines are 
isothermal SZ models with 8 C = 1/5 and AT(0) = 0.083 mK, 
(3 = 0.75; AT(0) = 0.050 mK, f3 = 0.5; and AT(0) = 0.12 mK, 
fi = 1.0, convolved as in (b). (d) Cross-correlation of WMAP 
94 GHz W data with 500 galaxy groups and clusters selected by 
friends-of-friends algorithm (fl = 8, m 3s 35) from the K s < 13.7 
2MASS XSC. In all cases the solid lines represent the average 
AT obtained over the area used. Monte-Carlo errors including 
the effects of clustering of clusters are shown. 



ters may have too low densities to produce a strong SZ effect. 
We conclude that there is a significant anti-correlation out 
to ~ 1 deg scales between the ACO R ^ 2 cluster catalogue 
and the WMAP W band data, confirming the tentative de- 
tection seen in the APM group and cluster catalogues, and 
that this is likely to be caused by the SZ effect. 

Figs. 2 (a-f) shows that the ACO R > 2 clusters 
also produce anti-correlations in all the other frequency 
bands, including the ILC. In terms of the predicted SZ fre- 
quency depend e nce acc ording to equations (11) and (13) of 
iRefreeier etafl fcOOOal) this is expected because relative to 
the W band the V, Q, Ka and K bands the SZ decrement 
should increase by the factors 1.16, 1.21, 1.25 and 1.25. The 
poorer resolution of the low frequency WMAP bands makes 
the detection of the decrements more unexpected but it is 
explained by the apparent persistence of the signal out to 
scales well in extent of the W band beam. The frequency 
dependence of the SZ signal appears statistically consistent 
with the SZ prediction (dashed lines) in all except the K 
band where the combination of poor resolution and residual 
Galactic contamination may be the cause of the poor fit. 
However, the errors are too large to discriminate between 
the SZ and CMB spectral indices. It is this spectral simi- 
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Figure 2. Cross-correlation of WMAP data in the W, V, Q, Ka, 
K and ILC bands with 606 |6| > 40 deg ACO clusters. The SZ 
decrement can be seen in each band. In all cases the solid line is 
the simple average AT over the area surveyed, and the vertical 
short-dashed line indicates the beam half-maximum in each band. 
The long-dashed line shows a S -1 / 3 fit to the W band, scaled for 
the SZ frequency dependence to the other bands, which can be 
used as a reference line for the frequency dependence. 



larity which makes it difficult to reject SZ contamination in 
the CMB maps and which explains our SZ detection in the 
'clean' ILC map. The APM clusters and groups also show 
anti-correlations in the V and Q bands but less so in the 
other bands, probably because of the lower resolution. 

To check if the cross-correlation results contain arte- 
facts characteristic of residual Galactic foreground contam- 
ination, we have cross-correl ated the ACO cluster s with the 
WMAP foreground maps of iBennett et al.l (I2003D and find 
no indication of any strong systematics. We also note that we 
have obtained results for the SZ decrements from the ACO, 
APM and 2MASS clusters from the 'foreground-cle aned' 
map of lTeemark. de Oliveira-Costa fc Hamilton] i2003T) that 
are entirely consistent with those found from the ILC map; 
the detection of the same anti-correlation in both of these 
'clean' maps further argues against the possibility that the 
anti-correlation is due to foreground systematics. 

As a further powerful check against systematics caused 
by foreground contamination, we have also applied the cross- 
correlation analysis after rotating the ACO R ^ 2 clusters 
around the Galactic Poles with respect to the WMAP W 
and V band data. Fig. 3 shows the W and V band cross- 
correlation signals (relative to the mean AT) integrated to 
the beam half-maximum (9 < 6. '3 for W, 9 < 9. '9 for V), 
9 < 60' and 9 < 500' measured at intervals of Al — 10 deg 
of Galactic longitude. The solid line shows the difference 
between the beam half-maximum and the 9 < 500' angular 
scales; this is to test for the persistence of the beam-size anti- 
correlation under the assumption that the 9 < 500' result 
is entirely due to systematic effects. The anti-correlation at 



Figure 3. (a) Cross-correlation of WMAP W band with the 606 
ACO clusters after rotation about the Galactic Poles by adding 
Al to the cluster's Galactic longitude. The open circles represent 
the average AT with 8 < 6/3 (the beam half-maximum in W), 
the crosses represent the average AT with 8 < 60' and the filled 
circles represent the average AT with 8 < 500'. In each case, 
the average AT over the area surveyed has been subtracted. The 
8 < 6. '3 and 6 < 60' results at zero rotation lie significantly 
below the others but the equivalent 8 < 500' result does not. The 
solid line represents the difference between the 8 < 6/3 and the 
8 < 500' results which shows the effect on the beam-size result if 
the anti-correlation observed at 6 < 500' is assumed due entirely 
to systematics. (b) As (a) for the V band. Here the open circles 
and solid line correspond to 8 < 9/9, the beam half-maximum in 
V. 



the beam half-maximum is clearly significant whether mea- 
sured relative to the overall mean (3.1cr for W, 3.6a for V) 
or the 9 < 500' results (2.8a for W, 3.1a for V). The anti- 
correlation at 9 < 60' also remains significant whether mea- 
sured relative to the overall mean (2.3a for W, 2.2a for V) 
or relative to the 9 < 500' results (2.0a for W, 1.7a for 
V). Although the rotated results appear more correlated at 
9 < 60' than at 9 < 6. '3, implying these latter significances 
be treated with caution, it remains the case that for the 
9 < 60' the zero rotation result shows the lowest AT. This 
is not the case for the 9 < 500' points where several of 
the rotated points show both lower and higher excursions 
around zero than the zero rotation result. Clearly there may 
be systematics which are beginning to dominate any real SZ 
signal at these very large scales. We conclude that the ro- 
tation experiments suggest that the observed temperature 
decrements around ACO clusters are significant from the 
beam-size out to 9 < 60' and these results are robust to 
the assumption that the cross-correlation is dominated by 
possible systematics on 9 < 500' scales. 

Fig. 1(d) shows the cross-correlation of the 500 2MASS 
groups and clusters limited at K s < 13.7 with the WMAP 
W band. Again the 2MASS data shows similar trends to 
the APM and ACO clusters with anti-correlation seen on 
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scales of the beam and extending at marginal significance 
out to scales of ~ f deg. Again the frequency dependence 
is consistent with an SZ signal (not shown). At the average 
2MASS group redshift, z w 0.06, 9 < 1 deg corresponds to 
r < 3k" 1 Mpc and thus here there is further evidence for 
possible extension of the SZ decrements to scales of order 
> l/i _1 Mpc. 



5 DISCUSSION 

We have found evidence for WMAP temperature decre- 
ments extending out to scales of 9 ~ f deg from galaxy 
groups and clusters, with hints that the decrements may 
extend to even larger scales. Previous cross-correlation of 
the ACO catalogue with the Rosat All-Sky Survey has 
shown diffuse X-ray emission from ~ f keV gas extend- 
ing out to ~ 2 deg from D = 5, R ^ 1 Abell Clus- 
ters, comp arable to the scale of the anti-co rrelation seen 
in Fig. lc iSoltan. Frevberg fc Hasineerll2002l) . Thus, based 
on this indication that the anti-correlation may be caused 
by the SZ effect originating from diffuse 'supercluster' gas, 
we make a first or der calculation of the overall Compton 
parameter, y. From IRefreeier et al.l (l2000af) we use the rela- 
tion ATsz/To = yj(x) where To is the CMB temperature, 
x = hv/kTo and j(x) is a spectral function which takes 
the value j(x) = — f.56 at 94 GHz. The simplest route is to 
go via the m ^ 7 APM groups and clusters because these 
have the biggest space density and so are the most rep- 
resentative of average sightlines. These have a sky density 
of 3.5 deg -2 with an average SZ decrement that extends to 
9 > 0.5 deg and so have a sky covering factor of approxi- 
mately unity. For 0.1 < 6 < 0.5 deg in Fig. lb, AT is rea- 
sonably flat so that effects due to the WMAP beam may be 
small and in this range ATsz = —3.0 ± 1.6 /j,K. Using the 
above relation with To = 2.726 K, t his converts into a value 
for y(z < 0.2) = 7±3.8x 10~ 7 . NowlRefreeier et al.N2000al) 
refer t o IScaramella. Cen fc Qstrikeri jl993h and jPersi^t^rt] 
( 1995) as suggesting that 40 per cent of y originates at 
2 < 0.2 in CDM models. On this assumption, we find 
y(z < oo) = 1.8±1.0 x 10~~ 6 . This compares to the 3a upper 
limit on the total integrated y parameter from the COBE- 
F1RAS measurement of th e spectral distortion of the CMB 
of ?/0 < oo) = 2.2 x 10" 5 jFixsen et al.lll996l) . The 3a up- 
per limit from cross-corr elating COBE DMR and FIRAS is 
y(z < oo) = 4.5 x 10" 6 jFixsen et alJll997f l. iBandav et alJ 
Jl996l) found a 3a upper limit of Sy(z < 0.2) < 1.5 x 10" 6 by 
cross-correlating COBE DMR with the ACO cluster cata- 
logue. Thus our result is not inconsistent with these previous 
observational upper limits. We also note that our estimate 
of y is 2-3 x higher than that p r edicte d in the SCDM model 
of IScaramella. Cen fc Ostrikeil (I1993T) and that it is similar 
to the value predicted in the ACDM model of iPersi et alJ 
(1995). 

We next fit isothermal models to the SZ decrements 
for the ACO R^=2 clusters. The mod els are taken from 
Equns. 15, 16 of lRefreeier et aD <|2000af) and convolved with 
Gaussians to represent the beams. We assume the value of 
(3 = 0.75 quoted for Coma and a value of 9 C — 1/5 which 
is the Coma value scaled approximately to the z — 0.15 av- 
erage redshift of the ACO sample. We find that 0.083 mK 
is the best fit for AT(0) compared to the 0.5 mK quoted 



for Coma. The results are shown in Figs. 1(b), (c). It can be 
seen that the data appears to show a more extended decre- 
ment than the model. The same model with (3 = 0.5 gives 
an improved fit at larger scales for AT(0 ) of 0.05 mK. 

T he cluster correlation function feahcall fc Soneiral 
1983) suggests the number of excess clusters at 20' < 9 < 
100' from an average cluster is « 1.3. Based on the average 
decrement seen at 9 < 20', we estimate this cluster excess 
will contribute AT = 0.5 uK, at 20' < 9 < 100', compared 
to the observed AT = 6 uK. Thus it does not seem possible 
for clustering of 'beam-sized' clusters to explain the extent 
of the decrement; an extended (9 > 1 deg) gas halo around 
individual clusters appears to be needed. 

We next calculate the gas mass associated with an ACO 
R ^ 2 cluster. Using equation (16) o flRefregier et al.l |2000a'l 
and assuming fcT = 5keV, r c = 0.2/i _1 Mpc, the W band fit 
with AT(0) = 0.083 mK gives a value for the central electron 
density, wp = 1.8h x 1 0~ 3 cm~ 3 . Integrating equn. (14) of 
IRefreeier etafl feOOOal) with (5 = 0.75 to r < 1.75ft -1 Mpc 
(« 13') then gives a gas mass of M « 3 x 10 13 /i~ 2 M Q , 
which is not unreasonable compared to the X-ray gas mass 
of M « 1 X 10 14 /t~ 2, 5 M(7) detected w ithin a similar radius 
in the Coma cluster <Lea et alJll97l . 

A model with /3 = 0.75, 9 C = 2. '3 and AT(0) = 
0.015 mK is shown to give a reasonable fit at small scales 
to the APM m ^ 7 sample in Fig. lc. 9 C = 2. '3 is the Coma 
value of r c = 0.2/i -1 Mpc scaled to z = 0.1. Again there is 
evidence that the profile is more extended than the model. 

We now proceed to calculate the average gas mass 
associated with an APM m ^ 7 group /cluster. Us- 
ing equation (16) of IRefreeier et alJ i2000al) and assum- 
ing r c = 0.2ft -1 Mpc, the W band fit with AT(0) = 
0.015 mK gives a value for the central electron den- 
sity, no = 1.6ft x 10~ 3 c m~ 3 (fcT /l keV). Integrating equn. 
(14) of IRefreeier eta!] feOOOal) with (3 = 0.75 to r < 
1.75/i _1 Mpc (w 20') then gives a gas mass of M w 3 x 
10 13 fr~ 2 (fcT/lkeV)M Q . Since the APM groups and clusters 
are more numerous than the ACO clusters, we choose to 
use them to make an estimate of Qo gaa . Assuming the spac e 
density of APM gr oups quoted bv lCroom fc Shanks! il999l) : 
iMvers et all J2003t) of 3 x 10" 4 /i 3 Mpc" 3 we find a value of 
fto 9as w 0.03/1" 1 (lkeV/fcT)(6» m£l;c /20') ' 75 . Assuming h = 
0.7, kT = 1 keV and 9 max = 20' gives Q. 3as w 0.04, close to 
the W MAP tt baryon = 0.044 ± 0.004 result (ISpereel et alJ 
2003). There may be evidence in Fig. 1(b) for the SZ decre- 
ments to extend to 9 m ax ~ 60'. In this case, the masses 
associated with APM groups rise by a further factor of w 2 
implying £lo 9as ~ 0.1, now a factor of w 2x higher than the 
standard value for fl baryon . 

Is there any possibility that SZ decrements in the 
WMAP 94 GHz data could contaminate the acoustic peaks 
measurement? The question is interesting because the SZ 
correlations appear to extend out to 9 ~ 1 deg, the loca- 
tion of the first acoustic peak in the power spectrum. We 
have therefore run some simple models where circular areas 
representing clusters with SZ decrements are distributed at 
random over simulated CMB fields, laid down from Gaus- 
sian random fields drawn from CMBFAST power spectra 
JZaldarriaea fc Seliakl200ol) . If the SZ clusters have the same 
3.5 deg -2 sky density as the APM m ^ 7 groups and we as- 
sume that they extend only to 9 < 0.5 deg with a decrement 
ATsz = — 3«K then the SZ power spectrum of clusters is 
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two orders of magnitude below the first peak at wavenumber 
I — 220. However, if we assume that these groups and clus- 
ters extend unevolved in their gas content past the z < 0.2 
APM limit out to z < 0.5, then the group + cluster sky den- 
sity will rise to « 50deg -2 . If we further assume values for 
the SZ extent and decrement which are at the high end of 
the ranges so far suggested, i.e. 6 < 1 deg and ATsz = 5 /iK, 
then an SZ peak appears in the power spectrum which is of 
order ~30 per cent of the amplitude of the first acoustic 
peak at / = 220. The consequences for the lower amplitude 
acoustic peaks at higher multipoles could be even more seri- 
ous; we note that the temperature power spectrum from the 
CBI experiment already shows an excess at I > 2000 which 
may be due to a strong SZ eff ect from individual clusters 
{Pearson et alJl2002l: Isilkll2002h . The above models might 
run up against the y(z < oo) = 4.5 x 10 -6 upper limit of 
iFixsen et al.l l|l997|) . However, it is clear that the question 
of how much SZ effects contaminate the primordial power 
spectrum is re-opened by the spatial extent of the SZ sig- 
nal found in our results. Higher resolution CMB data and 
deeper group and cluster catalogues are needed to constrain 
further the SZ contribution from z > 0.2 clusters. 



6 CONCLUSIONS 

We have found evidence for anti-correlation between WMAP 
data and the positions on the sky of galaxy clusters derived 
from the ACO, APM and 2MASS surveys. We interpret the 
signal as caused by the SZ effect, inverse Compton scat- 
tering of the CMB photons by hot gas in groups and clus- 
ters of galaxies. The signal may extend to ft! 1 deg scales 
around ACO clusters, implying they have extended gaseous 
haloes which may also constitute a diffuse gas component 
in superclusters. We estimate the mean Compton y param- 
eter associated with z < 0.2 APM groups and clusters as 
y(z < 0.2) = 7±3.8x 10~ 7 . This is not inconsistent with pre- 
vious upper limits and with expectations from CDM mod- 
els. We have also estimated the average R 2 ACO cluster 
gas mass assuming kT = 5keV and found this to be in 
reasonable agreement with X-ray observations of individual 
ACO clusters within the central radius, r < 1.75/i _1 Mpc. 
The average mass of APM groups+clusters is found to be 
M(r < 1.75/i _1 Mpc) w 3 x 10 13 /r 2 (l keV/fcT)M Q which 
gives D. gaa m 0.03ft" 1 (1 keV/fcT). For kT = IkeV and 
h = 0.7 this value for Qo 9as « 0.04 is close to the value 
of Q.Q haryon = 0.044 in the standard model dSpergel et alJ 
2003). But because the X-ray temperatures are likely to be 
lower than 1 keV and because the SZ decrements show evi- 
dence for gas haloes which extend beyond r = 1.75/i _1 Mpc 
then this value for S}o 3as could well prove to be a lower limit 
on the true value. 

We have briefly investigated how SZ contamination 
might affect the location and shape of the acoustic peaks 
in the WMAP temperature and find that although there is 
little effect from the temperature decrements found so far, 
if they persist out to z ~ 0.5 with the amplitude and ex- 
tent seen at z < 0.2, then even the first acoustic peak at 
the 1 deg scale could be significantly affected. Further cross- 
correlation analysis of deeper catalogues of groups and clus- 
ters will be needed to judge the seriousness of this potential 
SZ 'contamination'. 
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